Observation of dynamic localization in periodically-curved waveguide arrays 
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We report on a direct experimental observation of dynamic localization (DL) of light in sinusoidally- 
curved Lithium-Niobate waveguide arrays which provides the optical analog of DL for electrons in 
periodic potentials subjected to ac electric fields as originally proposed by Dunlap and Kenkre [D.H. 
Dunlap and V.M. Kenkre, Phys. Rev. B 34, 3625 (1986)]. The theoretical condition for DL in a 
sinusoidal field is experimentally demonstrated. 
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The quantum motion of an electron in a periodic po- 
tential subjected to an external field has provided since a 
long time a paradigmatic model to study fascinating and 
rather universal coherent dynamical phenomena. These 
include the long-predicted Bloch oscillations (BO) for dc 
fields 1], i.e. an oscillatory motion of the wave packet 
related to the existence of a Wannier- Stark ladder en- 
ergy spectrum, and the more recently-predicted dynamic 
localization (DL) for ac fields [2], in which a localized 
particle periodically returns to its initial state following 
the periodic change of the field. In recent years, BO have 
been experimentally observed in a wide variety of systems 
including semiconductor superlattices , atoms in accel- 
erated optical lattices 4], and optical waveguide arrays 
with a transverse refractive index gradient ^, J., 7] . DL 
is a phenomenon similar to BO which occurs when the 
electron is subjected to an ac field. The condition for 
DL, as originally predicted by Dunlap and Kenkre 2\ in 
the nearest-neighbor tight-binding (NNTB) approxima- 
tion and for a sinusoidal driving field E(t) = Fsm{ujt)^ 
is that Jo(r) = 0, where T = eaF/hoo and a is the lattice 
period. DL has been shown to be related to the collapse 
of the quasienergy minibands |^ , and the general condi- 
tions for DL beyond the NNTB approximation and for 
generalized ac fields have been identified 9]; DL under 
the action of both ac and dc fields has been also stud- 
ied and the influence of excitonic and many-body 
effects on DL in semiconductor superlattices has been 
considered (see, e.g. 11, 12j). Despite the large amount 
of theoretical work on DL, experimental evidences of DL 
are very few 0, and much less persuasive than those 
reported for BO. In semiconductor superlattices, the oc- 
currence of a variety of detrimental effects makes it 
hard an unambiguous demonstration of DL. The observa- 
tion of dynamical Bloch band suppression for cold sodium 
atoms in a frequency-modulated standing light wave has 
been related to DL Q, however a direct experimental 
evidence of DL in re al sp ace is still lacking. Recently, 
it has been suggested [l^ [l^ that optical waveguide ar- 
rays with a periodically-bent axis may provide an ideal 
laboratory system for an experimental realization of DL 



in optics, the local curvature of the waveguide providing 
the optical equivalent of an applied electric field |l^ IT3 • 
Since optical DL corresponds to a periodic self-imaging of 
light 16], it bears a close connection with self-collimation 
and diffraction cancellation effects that have been already 
demonstrated in photonic crystals (PC) and waveguide 
arrays 18^, 19, 20]. In these previous studies, diffrac- 
tionless propagation is achieved owing to the vanishing 
of diffraction at the inflection points of the isofrequency 
PC band surfaces |l^ |2^ or by alternating the sign of 
diffraction 18] using zig-zag waveguides (diffraction man- 
agement). 

In this Letter we provide the first direct experimental 
observation of DL of light in sinusoidally-curved waveg- 
uide arrays which exactly mimics the original Dunlap 
and Kenkre (DK) model 2]. The quantum-optical equiv- 
alence between light propagation in sinusoidally-curved 
waveguide arrays and the motion in a crystal of an elec- 
tron subjected to an ac field has been formally stated 
in Ref. 16]. The effective two-dimensional (2D) wave 
equation describing beam propagation in a periodically- 
curved array reads (T^ 
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where: z is the paraxial 
A = A/(27r) = l//c, V{x) 
the effective refractive index 
with period a [n{x + a) = 



propagation distance, 
^ Us — n(x), n{x) is 

profile of the array 
n(x)]. Us is the sub- 



strate refractive index, and 0:0(2;) describes the pe- 
riodic bending profile of the waveguide with period 
A A. By means of a Kramers-Henneberger trans- 
formation x' = X — xo{z), z' = and (\)(x' ,z') = 

i\){x' , z') exp -i{nslX)x^{z')x' - i{ns/2X) xg(^) 

(where the dot indicates the derivative with respect to 
z'), Eq.(l) is transformed into the Schrodinger equation 
for a particle of mass m = Us and charge in a periodic 
potential V{x') under the action of an ac field S{z'): 



iX 



dz' 



A^ 

2n. dx''^ 



+ V{x')(}) - q£{z')x'(j), (2) 



2 



where A plays the role of the Planck constant and the 
ac force is related to the axis bending by the equation 
q£{z') = —nsX{){z'). Note that, in the optical analogy, 
the temporal variable of the quantum problem is mapped 
into the spatial propagation coordinate z' ^ so that DL is 
simply observed as a return of the initial light intensity 
distribution during propagation. For single-mode waveg- 
uides, in the NNTB approximation and assuming that 
the lowest Bloch band of the array is excited, from Eq.(2) 
the following coupled-mode equations can be derived 

iCn = -A(Cn+l + Cn-l) - ^^£{z')Cn, (3) 

for the amplitudes of the field in the individual waveg- 
uides, where A > is the coupling constant. The phe- 
nomenon of DL 2] corresponds to periodic self-imaging 
at planes z = 0, A, 2A, 3A, .... In the NNTB approxi- 
mation, the condition for DL is L ( i£ exp[— Z7(^)] = 0, 
where 7(z) = {Ansa/ X)xo{z) 0,10. In particular, for 
a sinusoidally-bent array with period A and amplitude 
A, xq{z) = 74sin(27rz/A), the condition for DL is Q 
Jo(r) = 0, with 

AX ' 

The onset of DL condition Jo(r) = can be at best 
visualized by a direct measure of the impulse response 
[cn{0) = ^n,o] of a curved array As for the straight 

array the light spreads linearly with propagation distance 
z as |cn(2;)p = |Jn(2Az)p, yielding for the mean square 
number of excited waveguides the diffusive law (n^) = 
^^n^|cnp = 2A'^z\ for the sinusoidally-curved array, 
one has instead [2, l2l| 

{n^)=2A^[u\z)^v\z)], (5) 

where we have set u{z) = dr cos{Sor]{r)), v{z) = 
drsm{Sor]{r)), So = 47r^ans74/(AA^), and r]{z) = 
[1 — cos(27rz/A)]/[27r/A]. In particular, for a modula- 
tion period A smaller than the coupling length 1/A, one 
has \cn{z)\^ - I Jn(2AJo(r)z)p and (n^) 2A^J^{r)z\ 
Therefore, when the condition Jo{T) = is satisfied, an 
effective suppression of waveguide coupling is attained, 
which corresponds to a coherent destruction of tunnel- 
ing l2l|. For Jo(r) 7^ 0, the light diffraction pattern 
is analogous to that of a straight array with an effec- 
tive coupling coefficient Ag// = A|Jo(r)|. In our ex- 
periments, we could vary the parameter F [Eq.(4)] by 
either tuning the probing wavelength A or by measuring, 
for a given wavelength, the diffraction patterns of sev- 
eral arrays with different modulation period A. The first 
realized sample (Fig.l) consists of a set of one straight 
and two sinusoidally-curved arrays of waveguides with a 
graded index profile fabricated by the annealed proton 
exchange (APE) technique in 2:-cut congruent lithium 
niobate [2^. Channel width 7 /urn) and fabrication 
parameters have been chosen such that the individual 




FIG. 1: (color online) Left: Schematic of the first manufac- 
tured sample showing the three waveguide arrays: the straight 
array (1), the sinusoidally-curved array with short period (2), 
and the semi-cycle sinusoidal array (3). Right: Microscope 
image showing a particular of the straight and short-period 
curved arrays. 
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FIG. 2: Left: 2D output intensity pattern recorded on the 
IR camera for single waveguide excitation at A = 1610 nm of 
the straight array, and corresponding cross section intensity 
profile. The points show the waveguide intensity distribution 
|cn|^ = I Jn(2AL)|^ predicted by the coupled- mode equations 
(3) for A = 3 cm~^. Right: beam evolution (top view) pre- 
dicted by numerical simulations for the straight array. 



waveguide turns out to be single- mode, for TM polar- 
ization, in the whole spectral range from A = 1440 nm 
to A = 1610, with an upper cut-off frequency at around 
A ^ 1660 nm. Each array is made by a set of 80 identi- 
cal and equally-spaced L = 28-mm-long waveguides with 
separation a = 14 /im. Array (2) comprises 7 sinusoidal 
cycles (A = 4 mm and A c^^ 13 /im), with F ^ 2.405 at 
A 1610 nm, whereas array (3) comprises only a semi- 
cycle (A 56 mm and A 164 /im), with F 2.405 at 
A :^ 1440 nm. The waveguide coupling constant A has 
been measured as a function of the probing wavelength 
A by the analysis of the diffraction pattern obtained for 
the straight array (1) with a single waveguide excitation 
at the input plane (Fig. 2). Light coupling was accom- 
plished by focusing the circular diffraction-limited beam 
of a tunable semiconductor laser (Agilent Mod. 81600B) 
into the input face of one waveguide of the array to ob- 
tain a focused spot size of radius 3.5 /im. The light 
beam distribution at the exit of the array was imaged 
and recorded onto an IR Vidicon camera. The measured 
coupling constant A increases as the wavelength is in- 
creased and varies from A ~ 1.75 cm~^ at A = 1440 nm 
to A ~ 3 cm~^ at A = 1610 nm (the case shown in Fig. 2). 
The measured impulse response of the curved array (2) 
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FIG. 3: Left: 2D output intensity patterns recorded on the IR 
camera for a single waveguide excitation of the short-period 
curved array at different wavelengths, and corresponding cross 
section intensity profiles (solid curves). In the figures the 
dashed curves are the numerically-predicted intensity profiles. 
Right: corresponding beam evolution (top view) predicted by 
numerical simulations of Eq.(l). 
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FIG. 4: Same as Fig. 3, but for the semi-cycle curved array. 



for a few values of A is shown in Fig. 3, together with 
the beam evolution along the arrays as predicted by a nu- 
merical analysis of Eq.(l). Note that, as the wavelength 
A is tuned at the F = 2.405 value (top figure) a single 
spot is observed at the output plane, whereas far from 
the DL condition the number of excited waveguides in- 
creases. Note however that the change of F around ~ 2.4 
allowed by wavelength tuning is rather modest (about 
10%), which explains why the number of excited output 
waveguides is relatively small (~ 3 — 4) even at A = 1440 
nm. An important property of DL with a sinusoidal ac 
modulation is that wave refocusing is attained not only 
at multiples of the modulation cycle as shown for array 
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FIG. 5: Measured 2D output intensity patterns and cor- 
responding cross sections for an astigmatic Gaussian beam 
excitation of array (3) with (a) Wx 24.7 /xm, and (b) 
Wx ~ 37.4 fim. The dashed curves are the input Gaussian 
beam profiles in the x direction. 



(2) in Fig. 3, but also at multiples of the semi-cycle. This 
is clearly demonstrated both numerically and experimen- 
tally in Fig. 4, which shows the impulse response of array 

(3) . We note that, in the context of DL, the refocusing 
property in a semi-cycle was theoretically investigated 
for an ac square- wave driving field 15, 23] and related to 
periodic BO motion with alternating sign for the dc field 
(the so-called ac BO ^3)- However, the dynamics shown 
in Fig. 4 is not a BO motion with a dc field, which would 
strictly require circularly-curved waveguides 0, . 
We also checked the onset of DL by broad beam excita- 
tion illuminating at normal incidence the input facet of 
the array with an astigmatic Gaussian beam of varying 
width Wx in the in x direction. As an example. Fig. 5 
shows diffraction cancellation observed in array (3) with 
excitation at A = 1440 nm for two values of Wx] similar 
results are obtained for array (2). It should be noted 
that an analogous diffraction suppression was previously 
observed in Ref-ilS] using zig-zag arrays (see, in partic- 
ular. Figs. 3 and 4). In that case, diffraction cancella- 
tion arises from the alternation of positive and negative 
diffraction, obtained by periodic tilting of straight waveg- 
uide pieces. This is basically the same idea of "dispersion 
management" known for temporal pulse propagation in 
dispersive media. This scheme of diffraction cancella- 
tion, however, can not be regarded to as a realization of 
DL, i.e. wave localization induced by a periodic ac force, 
as intended within the original DK model (sinusoidal ac 
field) or its generalizations Q , including a square- wave ac 
field 15]. In fact, for zig-zag pieces of straight waveguides 
the curvature xq is zero and the £ field vanishes in the 
quantum- mechanical model [see Eq.(2)]. At the tilting 
planes of waveguide axis, however, the field wave vector 
(momentum) experiences an instantaneous change of its 
direction with respect to waveguide axis with a corre- 
sponding change of the diffraction sign, leading to a zero 
mean diffraction coefficient according to Eq.(4) of Ref. 

18]. In the quantum mechanical model (2), this would 
correspond to a quantum particle moving in a periodic 
potential whose momentum is abruptly and periodically 
kicked. 

To quantitatively check the diffusive law (5), we fabri- 
cated a second sample consisting of 7 sinusoidally-curved 
arrays with the same length and design parameters as 
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FIG. 6: Measured number of excited waveguides (n^) versus 
r (dots) for 7 waveguide arrays with increasing values of A, 
and corresponding theoretical behavior (solid curve). The 
recorded 2D output intensity profiles are also shown. 

those of the array (2) in Fig.l but with different and in- 
creasing values of A [from 2.8 mm to 14 mm], enabling 
to span a wide range of F values not achievable by simple 
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